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Continuing our analytic computation of the first-order self-force contribution to Detweiler’s red- 
shift variable we provide the exact expressions of the ninth and ninth-and-a-half post-Newtonian 
terms. 


I. INTRODUCTION 

The prospect of soon detecting the gravitational-wave 
signals emitted by coalescing compact binaries motivates 
a renewed study of the general relativistic two-body prob¬ 
lem. One of the useful lines of attack on this problem 
is the gravitational self-force program, which considers 
large-mass-ratio binary systems (mi <C m 2 ), and uses an 
expansion in powers of the mass ratio q = mi/m 2 1. 
Within this program, Detweiler [I| has emphasized the 
importance of focussing on the computation of gauge- 
invariant quantities, and he gave (for the case of circu¬ 
lar motions) the example of the function relating the 
redshift U* = dt/ds along the worldline of the small 
mass mi to the orbital frequency f2. To first order in 
q this gives rise to the gauge-invariant function Ul{y), 
where U*{y) = (1 — 3y)~^^‘^ q Ul{y) 0{q^) and where 
y = {Gm 2 ^/denotes a dimensionless parameter re¬ 
lated to the orbital frequency. Note that y can be consid¬ 
ered as measuring (in a gauge-invariant way) the dimen¬ 
sionless gravitational potential Gm 2 /(c^i?n), with Rq 
denoting the invariant radius canonically associated with 
n via Kepler’s law around the large mass: Gm 2 = CI^Rq. 
In the following, we shall denote the mathematical ar¬ 
gument of the first-order self-force function Ul{.) by u 
(evoking a gravitational potential) rather than y. [This 
notation is purely a matter of choice.] We henceforth also 
often set G = c = 1. 

Detweiler [I[ has shown that Ul{u) could be com¬ 
puted in terms of the first-order metric perturbation 
of a Schwarzschild metric of mass m 2 , say Sg^n, = 
g^j.„{x^-,mi,m 2 ) - g^^f^■, m 2 ) = qh^i,[x^) + 0{q^) as 

where 

( 1 . 2 ) 

Here denotes the Killing vector k = dt-\- and the 
superscript R denotes the regularized value of h^,j{x^) on 
the world line of the small mass mi. [When evaluating 
the first-order quantity /i^^. along a circular orbit of (co¬ 
ordinate) radius i?o it is enough to use the approximation 
Rq~ Rn = m 2 /u.] 


The beginning of the post-Newtonian (PN) expansion 
of the (first-order) self-force contribution Ul{u) was an¬ 
alytically derived in Ref. [Ij, namely : 

Ul{u) = -u-2u‘^ -5u^ - (1.3) 

Here, the first term {—u) is of Newtonian order, so that 
the second (—2 m^) and third (—5 m^) respectively repre¬ 
sent IPN and 2PN contributions. More generally, a term 
oc corresponds to the nPN level in Ul{u). 

The 3PN term was analytically derived (using full 
PN theory) by Blanchet et al. 0. In 2013, we Q 
showed how to analytically compute the 4PN term by 
a combined use of Regge-Wheeler-Zerilli (RWZ) formal¬ 
ism for the Schwarzschild perturbations together with 
the hypergeometric-expansion analytical solutions of the 
RWZ radial equation obtained by Mano, Suzuki and 
Takasugi Q (MST). We then progressively extended the 
analytical knowledge of the PN expansion of up to 
the 8.5 PN level [1,Q. [See the latter works for references 
to other related analytical studies.] 

Parallely to these analytical studies, Detweiler’s red¬ 
shift variable was numerically computed in Refs. P,|3)[3) 
and these numerical data were used to extract numerical 
estimates of several hi^er-order (then unknown) PN ex¬ 
pansion coefficients [l|,13 . A breakthrough in this extrac¬ 
tion of PN coefficients from numerical self-force calcula¬ 
tions was accomplished by Shah, Friedman and Whiting 
Q who numerically evaluated the MST hypergeometric- 
expansion of Ul{u) to one part in 10^^^ for orbital radii 
extending up to 10^°Gm2/c^ . This extremely high nu¬ 
merical accuracy on C/*(u) for extremely small values of 
the argument u allowed them to numerically extract PN 
coefficients up to the 10.5 PN level, and also to provide 
educated guesses for the exact analytical form of several 
high-order PN coefficients. 

We have shown inj^,Q that the results of Shah, Fried¬ 
man and Whiting Q agreed with our (fully) analytical 
ones up to the highest PN level we had then computed, 
namely the 8.5PN level. The aim of the present short 
note is to report on an extension of our analytical compu¬ 
tation to the 9.5 PN level (using the techniques explained 
in our previous papers), and on its comparison with the 
results of Shah et al. 
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II. NEW TERMS IN THE PN EXPANSION OF 
Ui{u) AT THE 9 AND 9.5 PN LEVELS 


Following the notation of Eq. (21) in Ref. Q, we write 
Ul{u) = 17i(m)|8.5pn + (ag - /Jg Inn + 79 ln^'u)'u^° 

+ (ag ,5 -/3g.5lnM)M^^/^ , (2.1) 

where C/((u)| 8 . 5 pn is known from Q. Here (ag, /3g, 7 g), 
and (ag. 5 , /^g.s) are the coefficients of, respectively, the 
9PN (m^°) and 9.5PN (u^^/^) terms, which we have now 
analytically derived. Our results for these coefficients 
are: 
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( 2 . 2 ) 


Our analytically derived results for yg and /Sg.s agree with 
the numerical-based analytical expressions previously ob¬ 
tained for these two particular coefficients by Shah et al. 
0. Concerning the other (newly analytically computed) 
coefficients, namely ag, /Jg, and Og.g we have indicated 
by boxes in the above equations the extent to which our 
results Mree with the numerical estimates given by Shah 
et al. Q. More precisely, the boxes above include one 
more digit than those given in Table I of Q. In all cases, 
the agreement is perfect modulo possible rounding effects 
on the last digit quoted in Q. 


HI. CONCLUDING REMARKS 


The analytic computation of the post-Newtonian ex¬ 
pansion of the first-order self-force contribution Ul{u) to 
Detweiler’s redshift function has been raised here 

to the nine and nine-and-a-half PN level, thereby provid¬ 
ing the exact analytical expressions of terms which were 
previously obtained only numerically by Shah et al. 

Let us finally note that, using the results of Refs. 
SIS, our results can be translated into the computa¬ 
tion of the nine and nine-and-a-half PN contributions to 
the linear-in-mass-ratio piece of the main radial poten¬ 
tial A{u] v) of the effective one-body formalism [l^ . 
Denoting, A{u-,iy) = 1 — 2u -I- z^a(u) -I- 0{iy‘^) (where 
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V = mim 2 /(mi + m- 2 )'^ = q/{l + qY), we have 

a{u) = a8.5PN(u) + (a^o + aYo Inw + a^o In^ u)u^° 

+K 0.5 + alo. 5 lnu)u 2 ^/^ (3.1) 

where a 8 , 5 PN(M) was given in Q, and where the newly 
derived 9PN and 9.5PN coefficients are: 
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The corresponding numerical values are (consistently with, but more accurately than in, Eqs. (27), (28) in [^) 


(3.2) 


aJjj = 4845.870557019444177347393421579822176656222929365... 
aYo = -8207.441915171961061495913671985378064622493890134... 
aYo = 178.777945237204496463755723014982274241533500792... 
aio.5 = -28324.307465213628065671194515396169336328722651715... 
aio.5 = 3275.813959067119914181314855451516851587807063565.... (3.3) 
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